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Abstract -We present a comprehensive analysis of effective temperatures based on fluctuation- 
dissipation relations in a model of an active fluid composed of self-propelled hard disks. We 
first investigate the relevance of effective temperatures in the dilute and moderately dense fluids. 
We find that a unique effective temperature does not in general characterize the non-equilibrium 
dynamics of the active fluid over this broad range of densities, because fluctuation-dissipation rela¬ 
tions yield a lengthscale-dependent effective temperature. By contrast, we find that the approach 
to a non-equilibrium glass transition at very large densities is accompanied by the emergence of 
a unique effective temperature shared by fluctuations at all lengthscales. This suggests that an 
effective thermal dynamics generically emerges at long times in very dense suspensions of active 
particles due to the collective freezing occurring at non-equilibrium glass transitions. 


Introduction. — Statistical mechanics provides a uni¬ 
fied theoretical description of systems at thermal equilib¬ 
rium in terms of the probability distribution over phase 
space, from which thermodynamic quantities such as tem¬ 
perature can be defined . A similar framework is lacking 
for out-of-equilibrium systems for which the definition of 
a temperature remains an open issue [^. Active matter 
formed by assemblies of living cell s [^ , bacteria [^, self- 
propelled colloids [5||^ or grains [^|10| , is a coherent class 
of non-equilibrium systems receiving increasing attention, 
because they raise fundamental issues and for potential ap¬ 
plications in soft matter and biophysics Eniii . A number 
of recent studies have addressed the question of whether 
“effective” thermodynamic concepts can be fruitfully ap¬ 
plied to describe the phase behaviour and microscopic 
dynamics of active matter. This question is natural be¬ 
cause if some mapping to an equilibrium situation exists, 
then the whole arsenal of equilibrium statistical mechan¬ 
ics becomes available for further theoretical treatment. 
In particular the definition of a non-equilibrium temper¬ 
ature 13 -18 , of an active pressure [8 19 20 , of activity- 
induced interactions [8 21 and non-equilibrium free ener¬ 
gies 22 23 have been investigated for self-propelled par¬ 
ticles. 

An effective temperature Teff defined as the parame¬ 
ter replacing the thermal bath temperature in fluctuation- 


dissipation relations was extensively studied in slowly re¬ 
laxing materials, such as spin and structural glasses 


27 


24- 


In these systems, the effective temperature becomes 

It can be used 


a meaningful thermodynamic concept 28 


to quantify heat flows and can be defined through a mi- 
crocanonical construction over a restricted region of the 
phase space. Physically, this approach is possible because 
there is a strong decoupling of timescales between ther¬ 
mal motion at short times, which essentially follows an 
equilibrium statistics, and an “effective” thermal dynam¬ 
ics at long times, which is an emerging collective property 
characterizing slow dynamical events leading to structural 
relaxation in driven and aging glasses. 

In active matter, the possibility to define a single¬ 
particle effective temperature in the dilute regime where 


particles do not interact has been discussed 14 16 17 29 


30 . The simple fluid state has been analysed mostly by 


numerical simulations determining fluctuation-dissipation 
relations for different observables in various models of 
self-propelled particles 13 15 18 . Finally, in the limit 


of very large densities where self-propelled particles may 
undergo a non-equilibrium glass transition, the driven 
glassy dynamics should resemble the one of slowly driven 
glasses and a eollective effective temperature was predicted 
to emerge from the analysis of a mean-field active glass 


model 31 , but this prediction has not been tested in a re- 
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alistic model in finite dimension. In this work, we focus on 
a model of an active fluid composed of self-propelled hard 
disks and analyse fluctuation-dissipation relations over a 
very broad range of densities encompassing the dilute, 
fluid and glassy regimes, in order to provide a compre¬ 
hensive picture of the relevance of effective temperatures 
in active matter. In particular, we find strong violations 
to the fluctuation-dissipation relations in all regimes, but 
conclude that effective temperatures appear most relevant 
in the slowly relaxing glassy regime. 


Numerical model. — We consider a two-dimensional 
system of N interacting self-propelled hard disks of di¬ 
ameter a, enclosed in a L x L square box with periodic 
boundary conditions. Self-propulsion is modelled by a 
non-Markovian stochastic drive, implemented by a kinetic 
Monte-Carlo (MC) rule that we describe below. The phase 
behaviour of the model and its microscopic dynamics were 
investigated before [8p2p3] . Here we extend these studies 
to analyse fluctuation-dissipation relations and effective 
temperatures over a large range of densities. 

A particle configuration is described by {ri{t),i = 
1 • • • N}, the set of particle positions at time t. The dy¬ 
namics proceeds as follows. At time t, a particle i is cho¬ 
sen at random and a small displacement di{t) of ampli¬ 
tude |di(f)| < So is proposed. Just as for equilibrium hard 
disks [^, the particle position is updated according to 
Viit + At) = Tiit) -I- <5i(<)Phard, where Phard = 1 if the 
move does not generate any overlap with a neighbouring 
disk, and Phard = 0 otherwise. 

Self-propulsion is introduced via a finite persistence 
time of the successive moves Siit). This non-Markovian 
dynamics breaks detailed balance. In practice, we choose 
5i{t + At) — Si{t) + r]^{t), where T]^{t) is a uniformly dis¬ 
tributed random shift of amplitude |i7i(t)| < Ji, chosen 
independently at each step. As a consequence, succes¬ 
sive displacements 6i decorrelate after a persistence time 
r = (i5o/<5i)^, and this dynamics generates persistent ran¬ 
dom walks for isolated particles. (See Ref. 


32 for more 


details about the model.) This model is appealing because 
it has only two control parameters. The persistence time 
r controls the self-propulsion of the particles, whereas the 
packing fraction ip = quantifies excluded volume ef¬ 

fects. In addition, the equilibrium hard disk model, where 
the physics is uniquely controlled by (/?, is restored in the 
limit T —>■ 0. The model does not include more complex 
features such as alignment rules or hydrodynamic inter¬ 
actions and serves as a minimal model to study the di¬ 
rect competition between glassiness and self-propulsion. 
A continuous-time version of the model for arbitrary par¬ 


ticle interactions has recently appeared 35 


We perform simulations using N = 10^ particles, vary¬ 
ing ip and T. We vary p from the dilute limit yi —>■ 0 up to 
p = 0.825, well beyond the equilibrium glass transition of 
hard disks p^^ « 0.80 36 , and r in the range r G [0,10^]. 


We fix for convenience do = cr/10. From now on, we set 
the Boltzmann constant ks = 1, we use a as the unit of 


length and one MC step represents N attempted moves. 
To suppress crystallisation in the glassy regime, we use 
a 35:65 binary mixture of hard disks with diameter ratio 
cti/ct 2 = 1.4 (cti then being the unit length) for systems 
above p = 0.69. For lower p, all systems are monodis- 
perse. Size polydispersity is introduced for convenience, 
but it does not influence the results presented below. 


Fluctuation-dissipation relations. — Using linear 
response theory at equilibrium, one can prove that the 
response of a system to an infinitesimal external pertur¬ 
bation is related to the spontaneous fluctuations of its con¬ 
jugate observable via the fluctuation-dissipation theorem 
(FDT) [^. Consider a system perturbed at time t = 0 
by an external field of constant amplitude /o coupled to 
an observable B. We define the cross-correlation func¬ 
tion between observables A{t) and B{t) as CAB{t,t') = 
{A(t)B{f))o — {A{t))o{B(t'))o, where (...)o denotes an 
equilibrium ensemble average. The conjugate linear re¬ 
sponse function RAB{t) is defined as 

{A{t)) - {A{t))o = fo f RAB{t,t')dt' + 0{fS), (1) 

Jo 

where (...) represents an average in the presence of the per¬ 
turbation. The time integral, XABit, t') = RAB(t, u)du, 
is the linear susceptibility, which is a more easily accessible 
quantity than Rab both in experiments and simulations. 
In terms of xab , the FDT becomes a simple linear relation 

XABitA') = Ji[CAB{t,t) - CABit.t')] ■ ( 2 ) 


Far from equilibrium, the FDT has no reason to be 
obeyed. However, it has been shown that in systems with 
slow dynamics, like spin glasses or supercooled liquids, the 
linear relation eq. ([^ is verified by replacing T by an effec¬ 


tive temperature Tf.s 1 ^ , 27 . Technically, one introduces 


the fluctuation-dissipation ratio XAB(t,t') through 24 


RAB(t,t') = ^ ^CAB{t,t')- (3) 


T 


If XAB{t,t') is slaved to the relaxation and becomes a 
function of the variable CABitA'), an effective temper¬ 
ature T^{C) = TIXab{C) can then be defined from 
the slope of a “fluctuation-dissipation (FD) plot” of xab 
vs. Cab parametrized by the time difference t — t'. 
Such plots have repeatedly appeared in the literature of 
glassy materials. The equilibrium FDT is recovered when 
Aab(C) = 1 and the slope of the parametric FD plot 
is then simply — I/T. The effective temperature defined 
above may in general depend on the chosen dynamic ob¬ 
servables A and B but in order to allow for a thermody¬ 
namic interpretation. Teg should be independent of this 
choice. This is of course the case in equilibrium and this 
hypothesis has been tested extensively in the glass litera¬ 
ture 


26 37 38 


We will investigate whether this important property also 
holds in our model. To this end, we first analyse the effect 
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of an external constant force /j = Cifoex applied from 
time t = 0 in the cc-direction, where = ±1 with equal 
probability. In our simulations, the perturbation is intro¬ 
duced via the acceptance probability of a MC update, 

^’acc(^z, /o) = min [l, Phard, (4) 

where Sxi = Si ■ and the equilibrium temperature is 
T = 1. For /o ^ 0, this introduces a bias in the x- 
direction, whereas the motion in the transverse y-direction 
is unaffected. We carefully checked in all our simulations 
that the value of /o we use is small enough to probe the 
linear regime so that deviations from the equilibrium FDT 
relations are uniquely due to the non-equilibrium nature of 
the self-propelled dynamics. The associated susceptibility 
is the displacement of the particles induced by /q: 


where x'^ = Xi + Sxi. A non-equilibrium effective temper¬ 
ature Teff 1 is then revealed if a straight line appears in 
the FD plot Xq vs. [1 — Fs{q, t)] for a given value of q. In 
the limit of small q, the long-time diffusion is probed and 
these observables reduce to the Stokes-Einstein relation. 


Effective temperatures in the active fluid. — In 

the dilute limit 0 , the only control parameter is the 
persistence time r. In this regime, particles have ballistic 
motion at short times, t < t, which crosses over to diffu¬ 
sive behaviour at long times, t > t, with D oc t 32 . The 


displacement xo induced by a constant force does not de¬ 
pend on the noise correlations, and hence is also indepen¬ 
dent of the persistence time. The mobility is then equal to 
the one of non-persistent Brownian particles. Thus, from 
the Stokes-Einstein relation we find 


1 ^ 

Xo{t) = hm^ — 

For these observables, the FDT in eq. ([^ reads: 


(5) 


XoW = — 


( 6 ) 


where Ax‘^{t) = N~^ ~ a^i(0))^)o is the mean- 

squared displacement. In the long-time regime Ax'^{t) ~ 
2Dt and Xo(^) ~ where D is the diffusion coeffi¬ 
cient and /i the mobility. At long times, the FDT yields 
the Stokes-Einstein relation, which is generalized out-of¬ 
equilibrium to 

Teff = P/m, (7) 


with an effective temperature that in principle de¬ 
pends on both T and q). 

In order to analyse carefully to what extend Teg depends 
on the choice of observables, we consider also Aq and Bq 


defined as 37 


N N 

^ 9(0 = ‘^'^^jCos[qxj{t)], ( 8 ) 

i=i j=i 


where Xj{t) = Tjit) ■ e^- With this choice, the associ¬ 
ated correlation function entering in eq. is the self¬ 
intermediate scattering function 

1 ^ 

and the corresponding linear susceptibility is 

Xq{t) = lim {Aq{t) - Aq{0))/fo . (10) 

Jo—nj 

This perturbation is physically meaningful, as its wavevec- 
tor dependence allows to probe the FDT at various length- 
scales. Numerically, we use 


PacciSxi, fo) = min 


g2ei/o[cos(qa:f)-cos(ga:i)] 


Thard, (H) 


Teff (r, -)► 0) = P/m oc T. 


( 12 ) 


This linear dependence on the persistence time obviously 
crosses over to a constant Teg — f 1 as t —)• 0 where equi¬ 
librium is recovered. Similar results were obtained before 


in a variety of active models 16,17 29 


Using the sedimentation of self-propelled particles, a 
meaningful Teg should also be accessible from the den¬ 
sity profiles (p{z) along the direction of the gravity field. 
This has recently been investigated experimentally in 
suspensions of Janus colloids and numerically using the 
present model. In the dilute limit, at the top of the sedi¬ 
ment, both the the Janus suspension and the self-propelled 
hard disk model were found to display an exponential den¬ 
sity profile, just as an equilibrium ideal gas at temperature 
Teff. Up to numerical accuracy, the effective temperature 
extracted from this static measurement is equal to the 
one obtained from the Stokes-Einstein relation. Thus, di¬ 
lute suspensions of self-propelled particles can be seen as 
a “hot” ideal gas with effective temperature Tgg oc r. The 
agreement between sedimentation and Stokes-Einstein re¬ 
lation is reasonable, as the gravity represents a biasing 
field acting at <7 = 0. Overall, this suggests that a single¬ 
particle effective temperature is useful to characterize di¬ 
lute active particles at large enough lengthscales, where 
details of the self-propulsion mechanism do not matter. 

We now turn to finite densities, where many-body ef¬ 


fects may modify the simple behaviour found in eq. ( 12 ). 
In Fig. I^a-c), we show the results of numerical simula¬ 
tions for D, fi and Teg for a broad range of t and tp values. 
When 79 —>• 0, the leading behaviour exposed above is con¬ 
firmed by the numerics. As discussed in detail in 32 , the 


diffusion coefficient D at finite density is a non-monotonic 
function of the persistence time, and a density-dependent 
optimum value for a persistence time T*{ip) appears where 
the diffusivity is maximal, see Fig. [^a). This maximum 
arises because self-propulsion simultaneously accelerates 
the dynamics of single particles, but it also leads to the for¬ 
mation of dynamic clusters where particles can be trapped. 
Physically, r* delimits the boundary between a homoge¬ 
nous fluid phase and a cluster phase characterised by the 
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Fig. 1: Diffusivity (a), mobility (b) and effective temperature 
(c) as a function of r for several packing fractions in the ac¬ 
tive fluid. In the dilute regime, D ~ r, /r ~ const, T^ff ~ ^/r 
(dotted lines). In the clustering regime, D ~ r“^, ~ 

Teff ~ -s/t (dashed lines in (a) and (b)). The effective temper¬ 
ature increases with r, but it decreases with ip. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

1 - Fsiq,t) 


Fig. 2: Parametric FD plots obtained in the active fluid for 
T — 9 and p = 0.747 for different values of the wavevector q 
characterizing the perturbation. Whereas the Stokes-Einstein 
value for Tefi is recovered in the g —>■ 0 limit, a different value 
of Teff characterizes each q, suggesting that the active fluid 
regime is not characterized by a single effective tempetature. 
Equilibrium FDT is shown as a full line. 


presence of fractal aggregates. Deep in the cluster phase, 
the diffusivity decreases as D ^ because it is con¬ 
trolled by the residence time of particles at the surface of 
the clusters [^. 


By contrast, our simulations indicate that the mobil¬ 
ity decreases monotonically with increasing r at constant 
density, see Fig. Sb). This simple behaviour arises be¬ 
cause self-propulsion does not enhance the mobility of 
isolated particles but it strongly affects its value in the 
activity-induced cluster phase. Asymptotically, we find 
that fj, ~ , suggesting that /i is dominated by the av¬ 

erage time spent by the particles inside the clusters. This 
behaviour can be understood by considering that parti¬ 
cles either diffuse freely between the clusters with a diffu¬ 
sion constant given by D{ip —>■ 0) oc r, or are kinetically 
trapped within clusters of average size n oc i/t 


32 


The 

time spent inside the clusters thus scales as t x n, so that 
the mobility scales as as observed numerically in 

the large-r limit. 

By taking the ratio of D and /r, we obtain the g —?► 0 
effective temperature shown in Fig. [^c). The dilute limit 
behaviour Teg ^ r is numerically recovered. At finite den¬ 
sity, we observe that Teg increases monotonically with r, 
and is thus not affected by the non-monotonic variation of 
the diffusion constant. However, the effect of the interac¬ 
tions at finite density is obvious as the effective tempera¬ 
ture grows much more slowly with r when > 0. In fact, 
the data shows that, for a fixed value of the persistence 
time, Teg decreases strongly with Lp, showing that inter¬ 
actions actually decrease the “agitation” of self-propelled 
particles. Indeed, in the strongly interacting regime, the 


combination of D ^ and p ~ suggests that 

Teg(T, (/?) ~ y/r. Therefore, the scaling with persistence 
time at finite density is qualitatively different from the 
linear behaviour found in the dilute regime. 


In a model of repulsive active dumbbells studied re¬ 
cently in the fluid regime, Tgg has a non-monotonic de¬ 
pendence on the activity that becomes more pronounced 
by increasing the density 18 . In our model, many-body 


effects alter the scaling of Teg with r in a simpler man¬ 
ner, monotonically reducing the effective “heating” of the 
system. 


We now analyse perturbations at finite lengthscales, to 
probe more deeply the thermodynamic meaning of the ef¬ 
fective temperature determined at large lengthscales from 
the Stokes-Einstein relation. In Fig. we show represen¬ 
tative parametric FD plots Xqi^) versus Fs{q,t) for our 
model at (fi = 0.747 with a fixed t = 9 and several values 
of the wavevector q. As expected, the correlation and re¬ 
sponse functions of the system deviate strongly from the 
equilibrium FDT. However, the FD plots all follow a con¬ 
vincing straight line. This seems to imply that for each 
wavevector, a well-defined value of the effective tempera¬ 
ture characterizes the long-time dynamics of the system. 
However, we also observe that the value of this effective 
temperature continuously depends on the chosen wavevec¬ 
tor, and Teg decreases when q is increased. This behaviour 
suggests that a unique effective temperature does not de¬ 
scribe the dynamics of our model in the active fluid regime. 

These observations extend to finite densities the results 
obtained analytically in the dilute regime of active Brow¬ 
nian particles, where it is similarly found that perturba¬ 
tions at finite lengthscales cannot be described in an ef- 
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Fig. 3: (a): Parametric FD plots for 5 = drr and t — 9 and 
different volume fractions ip approaching the non-equilibrium 
glass transition, (b) Effective temperatures extracted from (a) 
plotted as a function of the relaxation time . The horizontal 
dotted line corresponds to the estimated asymptotic value of 
Tefi = 7.2 as Tc — >■ 00. The dashed line corresponds to an 
empirical power law fit, T^s = 7.2 + 42t,7°'^. (c) FD plots for 
T = 9, ip = 0.819 and different values of q. All wavevectors 
yield a unique value of the effective temperature Tefr ~ 7.2 
(dashed line). 


fective thermodynamical framework 17 . This conclusion 


therefore applies to the entire fluid regime, where sen¬ 
sitively depends on the choice of conjugate observables. 
More broadly, these results echo recent studies showing 
that extensions of equilibrium concepts to active fluids is 
in many instances not straightforward 20p2p^29||39p0 


Collective effective temperature in the glassy 
regime. — We finally investigate the model at very high 
densities, approaching dynamical arrest from the fluid. 
Recent work supports the idea that dense assemblies of 
self-propelled particles can display a dramatic slow down 
of the dynamics sharing strong analogies with the glassy 
dynamics of particle systems in contact with a thermal 
bath [^|35,4l|- 4^. It was found in particular that activ¬ 


ity shifts the glass transition towards higher densities (or 
lower temperatures), but leaves the main collective proper¬ 
ties of the slow relaxation near the transition qualitatively 
unchanged [M|[4H[4^ . 

We compute the susceptibility Xq{i) and correla¬ 
tions Fg {q, t) associated to a modulated perturbation of 
wavevector q in dense suspensions of self-propelled disks. 
We focus first on the parametric FD plots for fixed q = 47r, 
r = 9 and several packing fractions from Lp = 0.747 (liquid 
regime) to tp = 0.821 (very close to the non-equilibrium 
glass transition ipc{T = 9) « 0.823). The results are dis¬ 
played in Fig. Similar results are found for other values 
of the persistence time. 

After a relatively short time corresponding to Fg {q, t « 
r) « 0.85 the FD plots for different packing fractions are 
well-described by straight lines with a slope that evolves 


slowly with ip. In the density regime investigated, Tgs de¬ 
creases from Teff « 15 to Teff ~ 7 over a density regime 
where the relaxation time increases by roughly seven or¬ 
ders of magnitude. The short-time period corresponding 
to the non-thermal part of the FD plot is a distinctive 
feature of active systems. In aging glassy systems, this 
crossover corresponds to the two-step decay of correlation 
functions with an intermediate plateau regime. Our sys¬ 
tem presents a two-step decay of Fg{q,t) as well 31 , but 


the plateau value of the correlation function plays no role 
in the FD plots. Instead, deviations from equilibrium are 
controlled by the timescale of the self-propulsion mecha¬ 
nism, an effective temperature appearing when t ^ t. As 
similar trend was noted in mean-field calculations 


41 


From the data in Fig.j^a) we extract the value of the ef¬ 
fective temperature and present its evolution with density 
in Fig.§b). To better appreciate the distance to the glass 
transition we simultaneously measure the structural relax¬ 
ation time Ta of the system defined as Fg{qmax,Ta) = 0.2, 
and represent T^s as a function of Tq,; qmax corresponds to 
the first peak of the static structure factor. These results 
show that Teff decays as dynamical arrest is approached, 
and we estimate that it approaches an asymptotic value 
Teff « 7.2 as Ta —>■ 00 , which clearly differs from the equi¬ 
librium value. Therefore we see that the persistence time 
T not only shifts the location of the glass transition to¬ 
wards large densities, it also fixes the value of the effective 
temperature at the glass transition. 

The crucial prediction of the mean-field theory of glassy 
dynamics is the observable-independence of the effective 
temperature. To test this, we compute Xq(t) and Fg(q,t) 
for several different wavevectors at fixed r = 9 and 
ip = 0.819. Remarkably, the results shown in Fig. are 
consistent with a unique effective temperature, as the FD 
plots for different q are almost superimposed. In particu¬ 
lar, the difference with the data at lower density in Fig.j^is 
striking. Numerically, the trend is that the g-dependence 
observed in the fluid becomes less and less pronounced 
as the density is increased towards the glassy phase, in 
perfect agreement with mean-field predictions. 


Conclusion. — The central result of the present work 
is the numerical confirmation in a finite dimensional model 
of self-propelled particles that a collective effective tem¬ 
perature emerges in active systems near a non-equilibrium 
glass transition. We have also shown that this physical pic¬ 
ture does not apply in the active fluid at small and mod¬ 
erate densities, where only perturbations at large enough 
lengthscales appear effectively “thermal”, whereas finite 
lengthscales perturbations are not well described by an 
effective thermodynamical treatment. 

Physically, the key aspect of active glassy systems that 
allows a meaningful definition of the effective temperature, 
is a clear separation of timescales between microscopic dy¬ 
namics at small timescales that is highly sensitive to de¬ 
tails of the self-propulsion and the structural relaxation 
taking place at much larger timescales. It is the long-time 
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sector for which the definition of an effective temperature 
appears meaningful. In addition, because the system un¬ 
dergoes a simultaneous kinetic freezing of all degrees of 
freedom near the glass transition, the effective tempera¬ 
ture becomes independent of the probed lengthscales and 
is thus essentially shared by all of them. 

* * * 

The research leading to these results has received fund¬ 
ing from the European Research Council under the Euro¬ 
pean Union’s Seventh Framework Programme (FP7/2007- 
2013) / ERC Grant agreement No 306845. 
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